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ABSTRACT 

We show that waves propagating in a transition disk can explain the short 
term temporal behavior of Cygnus X-l. In the transition disk model the 
spectrum is produced by saturated Comptonization within the inner region of 
the accretion disk where the temperature varies rapidly with radius. Recently, 
the spectrum from such a disk has been shown to fit the average broad 
band spectrum of this source better than that predicted by the soft-photon 
Comptonization model. Here, we consider a simple model where waves are 
propagating cylindrically symmetrically in the transition disk with a uniform 
propagation speed (c p ). We show that this model can qualitatively explain (a) 
the variation of the power spectral density (PSD) with energy, (b) the hard 
lags as a function of frequency and (c) the hard lags as a function of energy 
for various frequencies. Thus the transition disk model can explain the average 
spectrum and the short term temporal behavior of Cygnus X-l. 

Subject headings: accretion disks — black hole physics — radiation mechanisms: 
thermal — relativity — stars: Cygnus X-l 
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1. Introduction 

The X-ray emission of the well known Black Hole system Cygnus X-l shows two long 
term spectral states. The system spends most of its time in the hard state where the 
spectrum can be approximated as an power-law with an exponential cutoff around 100 
keV. The radiative mechanism generally invoked to explain this spectrum is unsaturated 
Comptonization of soft photons in a hot plasma. A model where the hot plasma is in the 
form of an corona on top of a cold accretion disk (Haardt & Maraschi 1993; Liang & Price 
1977) has been ruled out by spectral modeling of the spectrum ( Gierlinski et al. 1997; 
Nowak et al. 1999). The hot plasma could then be a hot inner accretion disk with the soft 
photons either arising from the cold outer disk (Shapiro, Lightman & Eardley 1976), from 
a pre-shock flow ( Chakrabarti 1997) or from internally produced Synchrotron photons as 
in the ADAF model (Narayan 1996). 

In the transition disk model (Misra & Melia 1996), the structure of the disk is obtained 
using the standard a-disk model of Shakura & Sunyaev (1973) without the assumption that 
the effective optical depth (r e //) is greater than unity. For high accretion rates (>~ 10 18 
g sec -1 ), r e ff turns out ( self-consistently) to be less than unity for the innermost regions. 
The disk is then forced to cool by bremsstrahlung self Comptonization which causes the 
temperature to be much hotter than that obtained by assuming black body emission. The 
temperature rapidly increase with decreasing radius in this region. Since the Compton 
y-parameter is always greater than unity ( although r e ff < 1) the local spectrum is a Wien 
peak. The sum over all Wien peak from the different radii (and corresponding different 
temperatures) gives rise to a power-law which matches the observed spectrum of Cygnus 
X-l ( Misra & Melia 1996). A similar model has been used to explain AGN spectra 
(Maraschi & Molendi 1990). Fitting broad-band (2 - 200 keV) spectrum of Cygnus X-l 
revealed that the transition disk model fits the data better than a simple Comptonization 



-4- 



one ( Misra et al. 1998; Misra, Chitnis & Melia 1998). In particular an ad hoc additional 
Wien peak {kT 50keV) feature required in the Comptonization fit to the spectrum ( 
Gierlinski et al. 1997) is not required for the transition disk model fit. 

Another way to differentiate between these theoretical models is by studying the 
timing analysis of this source. Recently RXTE observed Cygnus X-l and measured the 
power spectral density (PSD) for various energy bins, the time lags between energy bins as 
a function of frequency and energy and the coherence function (Nowak et al 1999a). These 
observations indicate that if the time lag is due to light travel time through the plasma, the 
size of the plasma has to be around 10 4 GM/c 2 (Kazanas, Hua & Titarchuk 1997) and the 
seed photon source has to be isotropic ( Nowak et al. 1999b). To avoid the energy problem 
associated with such a large plasma, Bottcher & Liang (1999), proposed a model where the 
variations are due to "blobs" of cool, dense matter which spiral inwards through the hot 
plasma. However, it is not clear whether these models agree even qualitatively with all the 
trends in the temporal behavior of Cygnus X-l as observed by RXTE. 

Nowak et al. 1999b showed that the time delays could be due to linear disturbances 
propagating through the plasma provided the speed of propagation is much smaller than c 
and photons of different energies are predominately produced in different regions. While 
such a scenario is not natural in the framework of a Comptonizing plasma, the transition 
disk meets the above requirements. In this letter, we show that waves propagating in a 
transition disk at least match qualitatively with most of the temporal observations. 

2. Waves in the transition disk 

A complete solution for the steady state disk structure in the framework of the 
transition disk model can be obtained using physical parameters which are the mass 
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accretion rate, the black hole mass and the the viscosity parameter a (Misra & Melia 1996). 
However, for simplicity and clarity Misra et al. 1997 used an empirical model to fit the 
average photon spectrum of Cygnus X-l. In the empirical model ( which differs from the 
transition model only in the parameterization form) the temperature profile is assumed to 
be of the form: 

kTocr n J m keV, (1) 

where r is the radius and J(r) = 1 — (rj/r) 1 / 2 . is the radius of the last stable orbit, taken 
to be Ti = 3R g , with R g = 2GM/c 2 . Here n and m are free parameters obtained by fitting 
the spectrum. The local radiative mechanism for such disks is saturated Comptonization 
(i.e. a Wien peak) since the local y parameter obtained is greater than unity (Misra et al. 
1998). Since the gravitational power per unit area F g oc R~ 3 J(R), the photon flux at earth 
can be written as (Misra et al. 1998) 

f E (E)= / ° F(E, r)dr oc E 2 / ° {kT)~ A exp(-E/kT) r~ 2 dr , (2) 

where E is the energy of the photon and r is the radius to which the transition disk extends. 
Misra et al. (1998) fitted this empirical model and obtained the best fit temperature profile 
T(r) with n —2.5, m w 1 and the maximum temperature to be 70 keV. This profile is 
shown in figure 1. 

The temporal behavior of the source is studied using the count rate as a function of 
time in a certain energy band. It is useful to define a response function (Nowak et al. 
1999b) for a particular energy band 

g(r)= F(E, r )A(E)dE (3) 

J Ei 

where E 2 and E\ are the limits of the energy band and A(E) is the energy dependent 
effective area of the detector. In this paper, we have used approximate values of A(E) for 
the RXTE detector. The normalized g(r) for some of the energy bands used in this paper 
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are shown in figure 1. The plot indicates where most of the photons contributing to the 
channel arise in the disk. We have assumed throughout this paper that the transition disk 
extends till 50r s . 

Following Nowak et al. (f999b), for simplicity, we shall consider waves propagating 
cylindrically symmetrically in the disk with a uniform propagation speed (c p ). Further 
only waves moving toward a sink in the origin are considered. The sink at the origin is 
represented as p s (r, i) = 5(r)p s (t) in the wave equation. In this formalism the Fourier 
transform of the observed light curve in an energy channel (s(t)) becomes ( Nowak et al. 
1999b, Morse & Feshbach 1953) 



where P s (f) is the Fourier transform of p s (t) and Gj is the Fourier transform of the relevant 
two-dimensional Green's function. This is solved to be ( Nowak et al. 1999b, Morse & 
Feshbach 1953) 



Several assumptions have been made in the above derivation. First, the propagation 
speed (c p ) is assumed to be a constant. Second, the sink is assumed to be located at the 
origin instead of and third the temperature profile (T(r) is assumed not to respond to the 
wave propagations. These assumptions highlight the limitation of the analysis to compare 
with observations. 




(4) 




(5) 



3. 



Results 



In figure 2 (a), we show the predicted time lags versus frequency between various 
energy bands compared with the observed values taken from Nowak et al. (1999a). Here, 
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c p = 3000r 9 / sees = 0.2c for a 1OM Black Hole which is the only free parameter in this 
analysis. Only for those values of frequency where the coherence function is near unity is 
considered (i.e / < 10 Hz). Considering the simplifying assumptions made in this analysis 
the results do qualitatively match the observations while there are discrepancies at higher 
frequencies. 

Crary et al. (1998) have used approximately 2000 days of data from the BATSE to 
estimate the time lag between 20-50 keV and 50 - 100 keV bands in the frequency range 
0.01 to 0.2 Hz. In figure 2 (b) we show the predicted time lag for the corresponding 
energy bands. The slope of the curve is similar to the observations in this frequency range. 
However, the observed time lag is 2 higher than the predicted one. This could mean that 
the wave propagation speed is generally slightly smaller than 0.2c. 

In figure 3, we show the predicted time lags versus energy for three different frequencies. 
The data are obtained from Nowak et al. (1999a). The logarithmic dependence for low 
energies is obtained for this model. Moreover, the predicted slope of the curves decreases 
with frequency matching the observed trend. Again at high frequency (ss 10) Hz there is 
significant deviation from the observations. 

The PSD cannot be predicted in the framework of this simple model since the Fourier 
transform of the sink term (P s (f) in eqn (4)) is unknown. Using the observed PSD we 
obtain P s (f) and present the result in figure 4 (a). Figure 4(a) also shows the predicted 
PSD„, if P s (f) was constant (i.e. white noise). It should be noted that in this model there 
is no exponential cutoff in PSD„, below 100 Hz. It would have been difficult to explain the 
observed PSD if there had been such a cutoff. Even though, the observed PSD cannot be 
directly inferred, the model does make predictions for the ratio of PSD in various energy 
bands since this is independent of P s (f)- In figure 4 (b) we show two calculated ratio as 
a function of frequency and compare it with the observed values. The ratios have been 
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renormalized to one at a frequency 1 Hz. For a wide range of frequency the ratios are close 
to unity which shows that the PSD is not sensitive to the energy band. One also finds that 
there is a remarkable correspondence of the calculated ratios with the observed values as a 
function of frequency and for different energy bands. Note that for frequencies greater than 
10 Hz, the present analysis is not adequate because of the lack of coherence. 

In the framework of this model a transition in the temporal behavior might take place 
when kr 12 ~ 1 where r 12 is the distance between the maxima of g(r) for the two energy 
bands under consideration. This would occur at a frequency f tr ~ c v j (27rr 12 ) ~ 477/r 12 ~ 15 
Hz for r\2 ~ 30r 9 (figure 1). It is tempting to associate f tr with the observed frequency 
where the coherence function drops below unity. At higher frequencies than f tr damping or 
non-linear effects may be dominant giving rise to this effect. 

In summary, waves propagating in a transition disk can qualitatively explain most of 
the temporal properties of Cygnus X-l. In the analysis only one additional parameter to 
the steady state disk model is used to explain the temporal behavior of the system. 

There are several observational features like the high frequency lag and the "shelves" 
in the lag versus frequency plot which is yet to be explained. The absence of such features 
maybe due to the simplifying assumptions (described in the end of section 2) used in this 
analysis. The transition disk models are generically radiation pressure dominated and 
hence are subject to secular and thermal instabilities (Misra & Melia 1996). However, 
these instabilities may be suppressed by non local effects like advection or convection. Fast 
moving radial waves as required in the above model may be produced in such a quasi-stable 
disk. 

The author would like to thank Ron Taam for useful discussions. 
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Fig. 1 — Temperature as a function of radius (solid line). The response function (g(r)) for 
various energy bins (dotted lines). From left to right the energy bins correspond to (a) 14.1 
- 45 keV, (b) 6.0 - 8.2 keV and (c) - 3.9 keV. 

Fig. 2 (a) — Time lags as a function of frequency between - 3.9 keV and (from top to 
bottom) (a) 14.1 - 45 keV, (b) 8.2 - 14.1 keV (solid line) ,(c) 6.0 - 8.2 keV and (d) 3.9 - 6.0 
keV. The dots are the observed time lags taken from Nowak et al. (1999a) between - 3.9 
keV and 8.2 - 14.1 keV. 

Fig. 2 (b) — Time lags as a function of frequency between 20 - 50 keV and 50 - 100 keV 
(solid line). The dashed line is the best fit to the observed values in the frequency range 
0.01 - 0.2 Hz ( Crary et al. 1998). 

Fig. 3 — Time lag as a function of energy for three different frequencies, circles: 0.3 Hz, 
squares: 1 Hz and triangles: 10 Hz. Points with error bars are data taken from Nowak et 
al. (1999a) for the corresponding frequencies. 

Fig. 4 — Top: (a) The observed PSD for energy bin 0-3.9 keV (solid line). The inferred sink 
function Ps(f) required to match the observed PSD (dashed line). The PSD if Ps(f) was a 
constant (PSD W ) which corresponds to white noise ( dotted line). Bottom (b): The ratio 
of the PSD. PSD so/t is for energy bin - 3.9 keV. PSD hard is for energy bins 14.1 - 45 keV 
and 6.0 - 8.2 keV. Solid lines are the computed values while dotted lines correspond to the 
observed ratios taken from Nowak et al. (1999a). 
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